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LONG-TERM  GOALS 

This  project  is  a  fundamental  study  of  the  basic  dynamical  mechanisms  involved,  and  the  consequent 
theoretical  modelling  approaches  needed,  in  the  generation  and  propagation  of  internal  solitary  waves 
across  the  continental  shelf  and  slope. 

OBJECTIVES 

There  are  two  principal  objectives.  The  first  is  to  develop  and  refine  amplitude  evolution  equations  of 
the  Korteweg-de  Vries  type  to  the  point  where  they  can  be  used  as  validated  models  for  the 
propagation  of  internal  solitary  waves.  The  second  is  to  undertake  a  major  re-examination  of  the 
generation  process,  using  a  combination  of  theoretical  and  numerical  analyses,  and  emphasising  the 
distinction  between  two-dimensional  and  three-dimensional  mechanisms. 

APPROACH 

Our  approach  is  to  develop  an  understanding  of  the  fundamental  dynamical  processes  involved  through 
a  combination  of  theoretical  analyses  and  numerical  simulations.  Our  research  group  comprises  post¬ 
doctoral  fellows,  research  students  and  international  collaborators  who  make  long-term  visits.  We 
maintain  contact  with  those  making  field  and  laboratory  observations,  with  the  aim  of  establishing  an 
ongoing  interactive  collaboration  on  data  interpretation,  model  development  and  validation. 

WORK  COMPLETED 

For  the  first  objective,  the  development  and  refinement  of  amplitude  evolution  equations  of  the 
Korteweg-de  Vries  type,  the  main  focus  to  this  point  has  been  on  understanding  the  role  of  cubic 
nonlinearity  vis-a-vis  that  of  quadratic  nonlinearity  in  several  contexts.  First,  we  have  obtained  a 
correct  asymptotic  derivation  of  the  coefficients  for  the  quadratic  and  cubic  nonlinear  terms  in  the 
extended  Korteweg-de  Vries  (eKdV)  equation  for  background  flows  which  allow  for  arbitrary  density 
and  current  stratification,  and  importantly,  allow  for  a  free  surface.  Second,  we  have  examined  the 
effect  of  various  frictional  processes  on  the  family  of  solitary  wave  solutions  of  this  eKdV  equation, 
using  primarily  asymptotic  techniques.  Third,  we  have  incorporated  topographic  forcing  into  this 
eKdV  equation,  and  examined  both  the  upstream  and  downstream  wavetrains  generated  by  transcritical 
flow  interaction  with  an  isolated  topographic  feature,  and  the  interaction  of  an  internal  solitary  wave 
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with  an  isolated  topographic  feature.  Fourth,  we  have  examined  the  solitons  generated  by  various 
initial  conditions  in  this  eKdV  model,  and  demonstrated  some  striking  differences  from  the  well- 
known  situation  for  the  KdV  model. 

For  the  second  objective,  that  is  the  utilisation  of  a  combination  of  theoretical  analyses  and  numerical 
simulations  to  re-examine  the  principal  generation  processes  of  internal  solitary  waves,  we  have 
focussed  first  on  the  development  of  a  suite  of  two-  and  -three-dimensional  numerical  codes  for  this 
purpose,  and  secondly  on  the  development  of  theoretical  models  of  the  Boussinesq-type,  which  can 
incorporate  arbitrary  stratification  and  topography  which  may  vary  in  both  horizontal  directions.  Our 
two-dimensional  spectral  codes  for  the  Euler  equations  have  been  extended  to  three  dimensions,  and 
are  now  being  tested  and  validated.  A  theoretical  model  suitable  for  describing  the  generation  of  the 
internal  tide,  which  allows  for  general  stratification  and  topography  has  been  developed,  and  is  now 
being  analysed.  Preliminary  results  show  that  this  model  can  reconcile  internal  ray  generation 
mechanisms  with  mode  scattering  mechanisms. 

RESULTS 

First,  we  have  examined  various  aspects  of  the  eKdV  model  for  the  propagation  of  internal  solitary 
waves.  For  the  present  discussion  this  is  given  by, 

A  +AA+\xAA+vA  A  +  \A  +  r(A)  +  / „(x)  =  0,  (1) 

t  A-  A-  A-  AAA 

where  A(x,t )  is  the  amplitude  of  a  representative  isopycnal  displacement,  t  is  a  time-like  variable 
describing  the  evolution  of  a  solitary  wave,  and  x  is  a  phase  variable  describing  the  shape  of  the 
solitary  wave.  The  coefficients  A,p,v,A  are  determined  by  the  linear  long- wave  modal  function,  which 
in  turn  depends  on  the  background  density  and  current  stratification.  The  expression  T(A)  is  a 
dissipative  term,  which  can  take  several  forms,  and  / (x)  is  a  term  representing  the  local  effects  of 
bottom  topography.  Explicit  expressions  for  the  key  coefficients  p,v,  X  ,  and  also  those  of  the  higher- 
order  linear  dispersive  term,  and  the  nonlinear  dispersive  terms  not  shown  in  (1),  have  now  been 
obtained  for  arbitrary  density  and  current  stratification,  and  for  a  free  surface.  One  of  the  most 
important  features  to  emerge  from  our  work  so  far  is  the  role  of  the  cubic  nonlinear  term  with 
coefficient  v  vis-a-vis  that  of  the  quadratic  nonlinear  term  with  coefficient  p  .  This  is  immediately 
evident  in  the  richer  structure  of  the  solitary  wave  solutions  supported  by  (1)  for  the  canonical  case 
when  the  coefficients  are  all  constants,  and  there  are  no  dissipative  or  forcing  terms,  when  the 
compared  with  the  corresponding  family  of  solitary  wave  solutions  of  the  KdV  equation  (i.e.  (1)  with 
only  quadratic  nonlinearity).  These  are  shown  in  Figure  1,  where,  without  any  loss  of  generality,  it  has 
been  assumed  that  p,X  are  both  positive.  Thus,  when  the  coefficient  v  of  the  cubic  nonlinear  term  is 
negative  (Figure  la)  we  see  that  the  solitary  waves  resemble  those  of  the  KdV  equation  for  small 
amplitudes,  but  for  large  amplitudes,  they  are  much  thicker  and  reach  a  limiting  amplitude  of-p/v  , 
known  as  the  “thick”  wave.  On  the  other  hand,  when  the  coefficient  v  of  the  cubic  nonlinear  term  is 
positive  (Figure  lb),  there  are  two  families  of  solitary  waves.  That  family  with  positive  polarity 
resembles  the  KdV  family,  but  that  with  negative  polarity  is  quite  different  and  in  particular  has  no 
small-amplitude  limit;  instead,  there  is  a  lower  bound  for  an  amplitude  of  -2p/v  ,  and  solutions  of  (1) 
with  lower  energy  are  represented  by  breathers,  that  is,  solutions  which  resemble  pulsating  solitary 
waves.  Given  that  observed  internal  solitary  waves  are  often  quite  large,  these  two  key  differences 
from  the  familiar  KdV  theory  are  likely  to  be  very  significant. 
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Figure  la  (top)  Solitary  wave  shape  for  negative  cubic  nonlinearity. 
Figure  lb  ( bottom )  Solitary  wave  shape  for  positive  cubic  nonlinearity. 


For  instance,  we  have  considered  the  effect  of  dissipation  on  a  solitary  wave,  using  various  forms  of 
the  dissipative  term  r(A)  representing  Newtonian  damping,  laminar  or  turbulent  boundary  layer 
damping,  or  damping  due  to  interior  turbulence.  For  a  negative  coefficient  v  (with  p,  A  both  positive) 
we  find  that  the  decay  of  a  “thick”  wave  can  lead  to  the  formation  of  secondary  wave  packets,  while 
for  a  positive  coefficient  v  ,  the  decay  of  a  wave  of  negative  polarity  leads  to  the  formation  of  breather 
states,  which  resemble  a  wave  packet.  Analogous  new  features  have  emerged  from  our  study  of  the 
interaction  between  internal  solitary  waves  and  topography,  modelled  in  this  case  by  (1)  with  the 
topography  represented  by  f(x) ,  and  with  no  dissipative  term.  In  some  circumstances  for  a  negative 
value  of  v  ,  the  interaction  can  cause  a  small-amplitude  KdV-like  solitary  wave  to  transform  into  a 
“thick”  wave,  while  for  a  positive  value  of  v  a  solitary  wave  of  negative  polarity  can  be  transformed 
into  a  breather.  Further,  in  the  trans-critical  regime,  we  find  that  when  the  cubic  nonlinear  tenn 
dominates,  the  upstream  solitary  waves  generated  by  flow  over  topography,  can  either  become  highly 
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irregular  when  the  coefficient  v  is  positive,  or  adopt  the  shape  of  an  isolated  monotonic  bore  when  this 
coefficient  v  is  negative.  It  is  pertinent  to  note  here  that  our  newly-derived  expression  for  the 
coefficients  in  (1)  indicates  that  for  realistic  oceanic  conditions,  the  coefficient  v  can  have  either  sign. 

Second,  motivated  by  recent  satellite  and  in  situ  observations  that  show  internal  solitary  waves 
emanating  from  point  sources  on  the  continental  slope,  often  submarine  canyons,  we  have  focussed  our 
present  research  work  on  the  development  of  theoretical  and  numerical  models  which  can  incorporate 
two-dimensional  variability  in  the  bottom  topography.  Our  three-dimensional  spectral  codes  for  the 
Euler  equations  are  now  developed,  and  are  being  tested  and  validated.  These  codes  have  the  potential 
to  provide  higher  spatial  resolution  than  most  current  codes,  and  importantly  are  non-hydrostatic,  so 
they  can  simulate  both  the  generation  process  and  the  subsequent  evolution  into  solitary  waves.  An 
asymptotic  theoretical  model  suitable  for  describing  the  generation  of  the  internal  tide,  which  allows 
for  general  stratification  and  topography  has  been  developed,  and  is  now  being  analysed.  The  basis  for 
the  development  of  this  model  is  a  novel  decomposition  into  vertical  modes  which  allows  for  arbitrary 
topographic  variation.  Importantly,  we  retain  the  free-surface  mode  in  the  model,  so  that  we  are  not 
only  able  to  describe  the  forcing  of  the  internal  tide  by  the  interaction  of  the  barotropic  tide  with 
topgraphy,  but  also  have  the  potential  to  determine  the  feedback  on  the  barotropic  tide.  In  its  simplest 
form,  when  it  is  linearised  and  made  hydrostatic,  we  are  using  it  to  describe  the  generation  of  the 
internal  tide,  for  both  one-dimensional  and  two-dimensional  topography.  Our  results  in  the  former  case 
show  that  this  model  can  reconcile  internal  ray  generation  mechanisms  with  mode  scattering 
mechanisms. 

IMPACT/APPLICATIONS 

We  anticipate  that  the  results  obtained  will  inform  the  scientific  community  about  the  structure  of 
internal  solitary  waves,  their  behaviour  under  such  environmental  impacts  as  friction  and  topography, 
and  the  processes  which  favour  their  generation. 
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